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1 Introduction 



1.1 Motivation 

The Cosmological Constant A [Q is presently at the epicentre of contemporary 
physics 1^ and plays a significant role in fields, such as: cosmology, astronomy. 
General Relativity, particle physics, and string theory. 

Experimental evidence has been growing suggesting that the Cosmological 
Constant has a small but non-zero value In particular, the results 

from high rcdshift Type la supernovae observations suggest an accelerating Uni- 
verse (positive Cosmological Constant) where the age of the Universe and Hub- 
ble's Constant H are approximately, ~14.5 billion years and 70 Kms~^ Mpc~^ , 
respectively. Other indirect experimental evidence (Turner Bahcall et al.Q) 
has put bounds on the energy densities ratio's of matter, dark energy and cur- 
vature, given by: Sl^ = STvGNPmatter/i^H'^) ~ 0-3, = §jj7 ^ 0.6 - 0.7, 
= —Kc^/{R{t)HY ^ 0, respectively. In these definitions Gm is Newton's 
constant, Pmatter dcnotcs the density of matter, K the curvature coefficient and 
R{t) the scale factor. 

We note here that The Cosmological Constant has been identified with a 
"dark" exotic form of energy that is smoothly distributed and which contributes 
60 -70% to the critical density of the Universe j|] (See Cohn and Carroll and 
Press |§ for a pedagogical review of the Cosmological Constant and Turner |^ 
for a review of Dark Matter and Dark Energy) . 

Recent data from the Boomerang experiment provide a firm evidence for 
the case of fiat Universe, fi^ ~ 0. 

It was suggested that the effects of the Cosmological Constant could be felt 
on a Galactic scale 1^ which would lead to the modification of Newton's inverse 
square law. This approach was used in order to try and explain the behaviour of 
the Galactic velocity rotation curves. The dynamical effects of A at megaparcec 
scale were investigated by Axenides et al [^j. 

Some of the theoretical justifications for a non-zero Cosmological Constant 
are discussed by Overduin et al [ p^ . 

The above arguments based on observations motivate a program of deriving 
exact solutions of General Relativity with a non-zero Cosmological Constant A 
and eventually apply these solutions to problems in cosmology and investigate 
dynamical effects of A at various astrophysical scale. 

The main objectives of this paper are: first to derive new exact solutions 
of General Relativity, for a non-zero Cosmological Constant, in the context 
of the inhomogeneous Szekeres-Szafron family of space-time line elements, and 
discuss at a preliminary level the effect of inhomogeneity in cosmology. Sec- 
ondly, to present in a systematic manner the cosmological implications of the 
homogeneous Friedmann, Robertson, Lemaitre, Walker (FRLW) limit of the 
above inhomogeneous metric. Details of the cosmological implications of the 
exact inhomogeneous solution will appear elsewhere. The material of this pa- 
per is organized as follows. In section 2 we derive a full exact solution of the 
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gravitational field equations with a non-zero Cosmological Constant in the inho- 
mogeneous class of Szekeres-Szafron models. The solution in a closed analytic 
form is expressed by the WeierstraB elliptic function. Our motivation for study- 
ing this particular class of inhomogeneous models is twofold. First the solution 
described above is new, and in general the inhomogenous cosmology is only a 
partially explored territory. Also the Szekeres-Szafrom metric is fairly general 
and consequently our solution of the resulting field equations represents the 
most general exact solution of General Relativity obtained to date. We note 
that the lack of general exact solutions of the equations of General Relativity 
has been an obstacle in obtaining the full realistic implications of the theory for 
cosmology. 

Secondly, a homogeneous Friedmann, Robertson, Lemaitre, Walker limit can 
be derived from this family and therefore a comparison between inhomogeneous 
and homogeneous cosmology can in principle be made. However, in this paper 
we do not attempt any such comparison. This will be a subject of a future 
publication. We discuss however, at a preliminary level where we expect the 
effect of inhomogeneity to appear. In section 3 we derive the FRLW limit and 
solve the Friedmann equation for the scale factor in terms of Weierstral3sche 
modular forms using the techniques of section 2. A formal expression for the 
Bubble's Constant is also derived. Section 4 provides the necessary mathemat- 
ical background for the use of the elliptic functions and elliptic curves in the 
solutions obtained. In section 5 we systematically investigate the cosmological 
implications of the exact solution of the Friedmann equation. The last section is 
used for our conclusions. We also include in appendices some of the formal cal- 
culations that were performed for this paper, and we emphasize that the same 
techniques are also usefuU for the solution of other important non-linear partial 
differential equations of mathematical physics. The present work is expected 
to be of interest to a wide audience, ranging from observational astronomers, 
relativitists, to theorists and applied mathematicians. We have therefore tried 
to point the reader to the appropriate literature for the necessary physics and 
mathematics background, and kept the mathematical formalism to the appro- 
priate level for making the presentation self-contained. 



1.2 Inhomogeneous versus Homogeneous Cosmologies 

As already mentioned in the introduction one of the reasons for choosing to study 
inhomogeneous cosmologies was to explore the relatively unknown territory of 
inhomogeneous exact solutions of Einstein's field equations. The Friedmann, 
Lemaitre, Robertson, Walker (FLRW) models are a realistic first approxima- 
tion to the geometry and physics of our Universe. By deriving and investigating 
exact solutions of the field equations of General Relativity with general inho- 
mogeneoous metric elements the theory will be prepared to meet any challenges 
posed by an inhomogeneous Universe. As it is emphasized by Krasinski in |p4| 
inhomogeneous solutions are very convenient in considering problems such as 
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the anisotropies in the cosmic microwave background radiation generated by 
density inhomogeneities in the Universe, evolutions of voids, the issue of singu- 
larities, possible effects of the expansion of the Universe on planetary orbits and 
others. 

Together the books, "Exact Solutions of Einstein's Field Equations" pl| , 
and " Inhomogeneous Cosmological Models" including references within, 

provide a comprehensive review of the presently known families of exact solu- 
tions to Einstein's Field equations. 

The Szekeres-Szafron family of solutions were chosen for study as they rep- 
resent a well-known inhomogeneous cosmological model, whose metric is given 
by: 

ds^ = dt^ ~ e2'3(-'?'^-'*)(da;2 + dy^) - e^^C-^?/.-*)^^^ (l) 

where the functions a (a;, y, z,t) , (3 {x, y, z, t) are to be determined by the field 
equations. 

The characteristics, properties and formal solutions for the energy densities, 
for this family of exact inhomogeneous solutions has been collated and docu- 
mented by Krasinski | ]T2[ . 

The Einstein equations with a non-zero Cosmological Constant, and a perfect 
fluid source energy momentum tensor T^^ have the form, 

Gp,u Rfiiy — -^gt^iyR = —{^9nu + (2) 

with 

T^u = [P + ~p)u^u^ ~ Pg^,^, Uf,^Sl, (3) 

here , P and p are the four-velocity, pressure and fluid energy density re- 
spectively. In this paper we study the model with Cosmological constant and 
pressureless matter (P = 0) in Eq.|[ To the best of our knowledge the exact 
solution of the non-linear partial differential field equations @ resulting from 
the metric Eq.(||), presented in this paper is new and represents the first solu- 
tion in a closed analytic form. Here we acknowledge the paper by Barrow and 
Stein-Schabes Q who were the first to derive a solution which was defined as 
a formal integral. In particular, as we shall see a{x, y, z, t), I3(x, y, z, t) and the 
energy density are expressed in terms of the Weierstrafi's modular forms. Conse- 
quently, our exact solutions for the inhomogeneous Szekeres-Szafron space-time 
element with a non-zero A lie on an elliptic curve. 

2 Exact Solution of the Gravitational Field Equa- 
tions for the Szekeres-Szafron Line Element 



4 



2.1 Approach, Definitions and Assumptions 

This section will outline the approach, definition and assumptions required to 
derive the general expression for a{x,y,z,t) and I3{x,y,z,t) from the Field 
Equations. The following sections will describe the determination of the Ein- 
stein's Field Equations in the Szafron form, and the derivation of the exact 
solution for the line element. 

We start by assuming that there exist co-ordinates which support the metric 
line element Eq.(|l|). 

The solutions are classified as either type /?' = or /3' 0, where /?' = 
dp/dz, depending on what initial trial solution is chosen. In this section the 
more general P' ^ type will be studied . 

2.2 Einstein and Szafron Form of the Field Equations 

Determining the Einstein field equations from the metric is straight forward but 
tedious. Expressions for the non-zero components of the Ricci tensor Rfj,^. and 
scalar R: Rqq, R22, R33, R are given in appendix A for completeness. 

The original solution for this metric was given by Szafron [|l4[ , his starting 
point was to introduce a pair of complex variables given by: 

^■.= x + iy, £,:=x~iy. (4) 

Following his approach we re-derived the field equations which are given 
below (due to an ordering difference of terms in the metric, our subscript three 
is equivalent to his one, i.e. Gfi = G33). The field equations are: 

G"o -2GI~G1^2 {a + 2'l3 + a^ + 2(3^) ^-k{p + 3P) = -(-2A + K~p), (5) 



-G§=/3'-/3'(d-/3) = 0, (6) 

Gl = 46-2/3/?^- + e-2"^'2 -2'$- 3/3^ = -(-kP) = -A, (7) 

G| = (-/?" + a' 13' - /3'2) 2e-^'^{a^j + a^af) - [a + P + + (3^ + a'p) = -{-kP) = -A, 

(8) 

Gi=ai+f3^-adP-a) = 0, (9) 
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\e^^Gl = + {a^f - 2p^a^ = 0, (11) 
where ' = d/dz, and ■ — d/dt, and 

:— da/d^ = -{da/dx — i da/dy), a-^ := da/dS, = -{da/dx + i da/dy). 

(12) 

Also we redefined the energy density and pressure by including the Cosmological 
Constant in the energy-momentum tensor: 



np :— Kp + A, 
-kP := A. (13) 

2.3 Exact Solution of the Szekeres-Szafron Line Element, 

Extending the work of Szekeres Szafron Q proved that /3j = 0, which 
allowed him to propose a trial solution for the field equations of the form: 

/3 = Log [$(i,z)] + z.(z,e,T)- (14) 

Since /3' ^ 0, Eq.(||) can be integrated to obtain: 

a = Log[/i(z, e z) + h{z, l)^{t, zy{z, ^, (15) 

It is then straightforward to show that the trial solution identically satisfies 
the rest of the equations (||) - with h = h{z), leading to two defining 
equations given by, 



e-" = A(z)ee + B{z)^ + B{z)^ + C{z), (16) 



where A{z)^ B{z),C{z),h{z), are arbitrary functions (A(z),C(z) are real func- 
tions, while B{z) is complex), P denotes the pressure and K{z) is an arbitrary 
real function determined by:[J 

A{z)C{z) - B{z)B{z) ^ \[h-\z) + K{z)]. (18) 



^Szafron used h{z) = 1 
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Equation can be integrated to give : 



1 K f 

$2 = -Kiz) + 2Af (z)/* - -- y P^dt. (19) 

The solution for <i> gives us a {x, y, z,t) ,13 {x, y, z, t) and the general expres- 
sion for the metric. In this paper by putting A = — kP, see Eg. (|l3|) , and per- 
forming the integration in the above equation, we get the following non-linear 
partial differential equation: 



$2 = -K{z) + 2Af(z)/$ + -A<I>2. 

3 



(20) 



Equation (gO|) resembles the equation of an elliptic curve (see section 4) and 
therefore can be integrated parametrically by the use of WeierstraBelliptic func- 
tions |^(see appendix B). The solution for $ is: 



<^{t,z) 



M{z)/2 



p{u + e) - p(uo) 
where the Weierstrafi invariants are given by 



where p(wo) '■= ^ 



K{z) 
12 ' 



(21) 



and the time t is obtained from t = J ^du with the result: 



(22) 



t+.nz) = - 



(t(u + e~ Vq) ^ , ^ ^ 



M(z) 



(23) 



or alternatively 



* + /(^) = /f 

Here p[z\uj,uj') , C(z) = ({z\tJ,Lj'), and (t(z) = cr(z|cj, cj') denote the Weier- 
strafi' family of functions, Weierstrafi, and Weierstrafi zeta and sigma functions 
respectively. Also e is a constant of integration and f{z) is an arbitrary function 
of z. The definitions and properties of these functions will be discussed in detail 
in section 4. The matter/energy density is obtained by the first field equation 
. We also calculated p'(wo) from the defining equation for the elliptic curve 

{p'M? = 4p3(„o) _ g^pivo) - .93 (25) 

■^For other examples of elliptic integrals see |16[ 



log 



a{u + e — vq) 
g{u + e + Wo) 



+ 2{u + e) C(fo) 



(24) 
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substituting p(wo) = ^-^^ in ( p5| ) we can show that that (p'(wo))^ = ti^^^^C^); 
and @ can be rewritten as: - y^^^^^^l^^^. 

Substituting our expressions for a (x, y, z, t) and /3 (x, y, z, t), given by: 



= $(t,z)e^(^^^'^\ (26) 



,2: 



dLog{^{t,z)) ^ dv 
dz dz 



(27) 



into the metric, Eq.(h|) the final solution becomes: 



{p{u + e) - pivo))'^ {p{u + €) - p{vo)Y dz d 



z 



(p(?x + e)-p(i;o))2 ^ ' 

The above solution represents the most general exact solution of inhomoge- 
neous relativistic cosmology obtained to date. 

2.3.1 Solution for the Matter Density in terms of the Cosmological 
Constant 

In this subsection we derive the expression of the energy density p that con- 
stitutes a solution of the Einstein field equations within the Szekeres-Szafron 
family of space-time elements. 

The mass density p, is obtained from the first field equation (see Eq. (||)), 

Gl~2G\-Gl = -K{p + iP), (29) 

Substituting the solution of a, (3 in Eq.(||) (see Appendix C for details) the 
general expression for the energy density p is given by: 

(30) 

Redefining the energy density as in Eq.(|l^), the final expression for the 
matter density is given by: 

^ 2M., + 6Mi/,, 
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The Friedmann, Lemaitre, Robertson, Walker (FLRW) limit of the above 
expression will be discussed in the following section. Here we will only state 
that in this limit Kp cx -pr^jp-, where R{t) is the standard scale factor in the 
FLRW line-element. 

We now discuss at a preliminary level where we expect the effect of inho- 
mogeneity to be important. One of the areas where the exact inhomogeneous 
solution of the field equations of General Relativity derived in this paper may 
play a vital role is in the field of gravitational lensing. In particular as is 
discussed in some length in ]l7[ | light propagation through an inhomogeneous 
universe differs from that through a homogeneous one. Light propagating in an 
inhomogeneous universe is deflected due to gravitational fields, thus limiting the 
accuracy with which positions of distant sources can be determined. The basic 
equations governing the change of cross-sections and shapes of ray bundles had 
been discussed by R.K. Sachs There it was shown how light propagating 
through a region of space where the density is lower than the average density of 
the universe diverge relative to corresponding light bundles in a homogeneous 
universe of the same mean matter density , |^ , . R. Kantowski consid- 
ered light propagation in the 'Swiss-Cheese' cosmological model ||2^. Owing 
to the lack of exact solutions of the equations of General Relativity that can 
describe an inhomogeneous universe, it was assumed that it can be described on 
average as a FLRW universe whose metric is locally perturbed by matter inho- 
mogeneities ( although as discussed in how this works is not obvious at all, 
due to the non-linearity of the Einstein equations). For the limiting case, where 
a light bundle propagates far away from all matter inhomogeneities C.C. Dyer 
and R.C. Roeder obtained a modified redshift-distance relation ||r^. The exact 
inhomogeneous solution obtained in this paper for the metric and the matter 
energy density may be used for realistic investigations of light propagation in 
inhomogeneous space-times. 

Another interesting application of the exact solution derived in this paper 
is to address the problem of Galaxy formation. Lemaitre predicted that 
an inhomogeneous Universe, with a non-zero Cosmological Constant, should 
contain collapsing regions spread over an expanding background (FLWR). In 
his paper he suggested that the expansion of the outer regions were caused by 
the cosmological constant repulsion while "local" instabilities could provide the 
mechanism for the formation of galaxies [^7| . 

The above physical applications of the exact solution to the inhomogeneous 
cosmology with the cosmological constant will be the subject of a separate future 
publication. 

3 Homogeneous Cosmology 

3.1 Recovery of the Friedmann Equation 

In this section we will show how the Friedmann equation (FLRW limit) can be 
recovered from the Szekeres-Szafron metric by judicious choice of the various 
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functions, parameters and constants (See Krasinski M). 

The starting point is the metric Eq.(^, where e^^^'^'*^ ^QOL(x,y,z,t) ^^^.g given 
by Eqns.(p6|),(|2^) respectively and e^^ by: 

e-" = A{z){x^ + + 2Bi{z)x + 2B2{z)y + C{z). (32) 

Here A[z),Bi{z),B2{z),C{z) and h{z) are arbitrary functions while K{z) is 
determined by: 

AC-Bl-Bl = \[h-\z) + K{z)\. (33) 
Now by choosing: 

$(t, z) = /(z)i?(t), if(z) = Kof{z), M{z) = Moz^, (34) 

the FLRW limit |l2j in the Goode & Weinwright form can be derived, 
namely: 

ds"^ = dt^ - R^{t)[e^''{dx^ + dy^) + W^fv^^Jz\ (35) 

here, 

= {e-Kf)-\ (36) 
e"^ = /(z)[a(z)(a;2+j/2)_^2&(z)a; + 2c(z)y + d(z)]-\ (37) 

e, iiT are arbitrary constants and R{t), f{z),a{z),b{z),c{z),d{z) arc arbitrary 
functions subject to the constraint ad — — = e/4. 

If in addition we choose, Bi = B2 = 0,C — 1, A — j, and f{z) = z, then the 
FLRW limit is recovered, leading to the standard Friedmann equation, given by: 



- -^Pr^aUerR' + "-f ~ Kc\ (38) 

and the matter energy density eq. (|3l|) , which is a solution of the gravitational 
field equations, in this limit becomes: 



where Mq is a Friedmann mass integral ||2J]. In Eq. (|3^) K and A denote 
the curvature coefficient, and Cosmological Constant, respectively. 



Now by choosing suitable parameters, namely. 
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2m = — - — X constant, 



K' = Kc^, 



A' = Ac2 



(40) 



equation ( |3q ) can be rewritten in the standard Friedmann form (Eq.(£0|)), 

2M A' 



whose solution is given by: 



(41) 



R{t) 



M/2 



p{u + e,T) + K'/ 12' 



t 



-c\Ij 



Log 



a{u + € - vo) 
_a{u + e + vo) 



+ 2u CM 



and the cubic invariants are given by 

92 



1 



12 216 12 



' Ji/r2 



and p(wo) = - ^ 



12 ■ 



(42) 



(43) 



Alternatively, the solution for the scale factor and the time in terms of Jacobi 
theta functions are given by: 



1 

4^ 



^„+i(0) ^1(7) 



12 



Log 



(7 + 7o) 



(7o 



(44) 
(45) 



where 7 = ' 7o ~ ^ are the three roots of the given cubic equation. 

3.1.1 Theoretical Expression for the Hubble Constant 

The theoretical expression for Bubble's Constant is given below in terms of 
Weierstrafi, or alternatively, Jacobi theta functions. 

The solution of the Friedmann equation for the scale factor is given in para- 
metric form by eq.(p^, and the derivative of the scale factor with respect to 
time i? := ^ is given by : 



p'ju + e) 
p{u + e) + % ' 



(46) 
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where p'(z) = One can obtain ( ^6| ) by directly differentiating eq.(^2|) 

with respect to time. Alternatively by substituting ( ^ ) in cq.(p8[) and using 
the equation for the elliptic curve: 

{p'{u + e)f = Ap\u + e) - 52 p{u + e) - 53, (47) 

one rederives (|4^). This procedure provides a consistency check for our calcu- 
lations. 

The theoretical expression for Hubble's parameter for a homogeneous cos- 
mology is given by: 



We note that Hubble's Constant is proportional to p'{u). Equation ( ^ ) is 
remarkable since it connects the expansion rate of the Universe, a dynamical 
parameter, with the cosmological parameters and constants of nature in a highly 
non-linear way. 

Alternatively, Hubble's Constant can be expressed in terms of Jacobi theta 
functions, as follows: 

2 1 / 02 (7)^3(7)^4(7)^^(0) A .^g. 

^ ' M 4tj3 V O2 (0) 6z (0) 6i (0) (7) / ■ 



Also the deceleration parameter q = is given by the following expres- 

sion: 



p"iu + €)[p{u + e) + ^] 

1 = ._,9/.. , ,^ 1' 



P'2(m + e) 

[6p^-y2Mu + e) + ^] 
p'2 (?i + e) 



- 1. (50) 



4 Roots of the cubic in terms of the cosmologi- 
cal parameters and regions of the fundamen- 
tal period parallelogram that correspond to 
physical solutions 

In this section we shall discuss at which regions of the complex plane the 
Weierstrai3 function, p{z), and its derivative p', are real. We restrict atten- 
tion in these regions of the complex plane since we require that the scale factor 
is positive and real R > 0), and that the time and Hubble's Constant are real 
physical quantities. Due to the above necessary conditions it is important to 
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discuss the relevant properties of the elliptic functions appearing in our solu- 
tions. For further formulas and mathematical background the reader is referred 
to the book by Abramowitz |^ and the book by Silverman and Tate . 

The determination of the cosmological quantities such as the scale factor R, 
involves finding the roots of the cubic polynonial, —g2Z — 53, of the elliptic 
curve for various values of the cosmological parameters A, M and K. 

The three roots of the cubic, et , i = 1,2,3 can be determined by using 
the algorithm developed by Tartaglia and Gardano and their expressions in 
terms of A, M, K are given by : 



1/3 

1/3 



1 r 1 1 

ei = — |^^"18Af^^ + 6V-if^AM2 + 9M4A2| 

+ ^ - ISAf^A - &^J-K^kNP + 9Af4A2|^ 

o r / 1 1/3 

62 = Y^|ii:3-18A'/2A + 6V-if3AM2 + 9M4A2| 

(P f . 1/3 

+ 1^'^ - 18i\/^A - Q^J-K^KM^^ + 9Af4A2| , 

1/3 

+ |i^3 _ isM^A - 6V-^3AM2 + 9M4A2 J , (51) 



where p 



„27ri/3 



In general there are three types of solution depending on whether the dis- 
criminant is either zero, positive or negative. The solutions are outlined below: 

1. A = 0, the roots are all real and at least two of them coincide. The 
curve is no longer an elliptic curve, since it now has a singular point. The 
solutions are not given by elliptic functions and do not have modular properties. 



2. A < 0, the roots eiand 63 form a complex conjugative pair with 62 real. 
These solutions are given by elliptic functions and have modular properties. 

3. A > 0, all the roots are real and the solution is expressed in terms of 
elliptic functions which have modular properties. 

The lattice L of periods 2lo, 2io' associated to the given cubic equation (that 
corresponds to a choice of the parameters A, M, K) is determined by integrating 
certain Abelian differentials — described below. 

y 

Consider the case A < 0. In this case two of the roots 61,63 are conju- 
gate complex and one 62 is real. There is a complex conjugate pair of semi 
fundamental periods, lo and lo' , which are given by: 
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oo 



dz 




oo 



dz 



^4^3-322-53 



LU = 



(52) 



For A > 0, all the three roots 61^62, 63 of 4^3 — g^z — 53, are real and if the 
Ci are ordered so that ei > 62 > 63 we can choose the periods as: 



In this case lo and lo' are real and totally imaginary respectively. The period 
ratio r is defined by r = cj'/o;. 

The three roots of the cubic can be expressed in terms of the Jacobi theta 
functions p5[ |, and therefore have non-trivial transformation properties under 
the modular group SL{2,Z). For instance under the modular transformation 
T — ^, ei — > 63, 63 — > ei and 62 is invariant. The roots are modular invariant 
under subgroups of SL{2,Z). In fact the factor group PSL{2, Z)/T2 where 
r2 is the modular group of level 2 is isomorphic to S3, the permutation group 
with six elements acting on the letters 61,62,63. One can construct modular 
invariant quantities using the above transformation properties. The combination 
of cosmological parameters given by the relation (62 — e3)(ei — 63)""'^ is the 
absolute invariant of r2 . 

4.0.2 The Fundamental Period Parallelogram for A < 

The Fundamental Period Parallelogram (FPP) for a typical solution, A < 0, is 
shown in figure |. The fi gure also shows the position of the real root 62 , and the 
two complex conjugatives ei and 63. In order to obtain physically meaningful 
(i.e. real) values for the scale factor R, we restrict ourselves to a set of solutions 
for the Weierstrafifunction and its derivative that lie along the diagonals of the 
FPP. Mathematically this can be described by: 

The Weierstral3 function p (z) = p (u + e) , z = u + e, where e is a constant 
of integration, 

is real when the argument z — Re (z) + ilm (z) — m{Lu + lu') + n (a> — to') , 
with, m,n € R. 

Note that the WeierstraB function and its derivative are real for the diagonal 
that goes left to right, and real and totally imaginary, respectively, for the right 
to left diagonal. As z varies along diagonals of period parallelograms from to 
2{uj + uj') p{z) decreases from -\-oo to 62 to 00. 

We also note that the zeros of the derivative of the Weierstraf? function 
p'{x) = 2yJ{p ~ 61) (p — 62) (p — 63) are at the point lo , lo' , uj + lu' of the FPP in 
the complex plane. From the Euclidean uniformization f : u ~> (piu), p'{u)), u 




(53) 
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a complex number, of the Weierstrafi elliptic curve these points correspond to 
the three points of the eUiptic curve of order two n. 



4.0.3 The Fundamental Period Parallelogram for A > 



The Fundamental Period Parallelogram (FPP) for, A > 0, is shown in figure 
(|). There are three real roots 61,62,63 e R, and uj,uj' are real and totally 
imaginary, respectively. The real values of the Weierstrafi function lie on the 
vertical and horizonal lines defined by: 

z = u + e, p{z = a + ib) E R 

for a = or w and h — Q lo' . 



The Weierstrafi derivative is real on the horizontal lines and totally imaginary 
on the vertical ones. More specifically, in the interesting case of 32 — "^2" 7^ 
0, (73 = 0, the resulting lattice is a square lattice with uj' — ico, i.e., the lattice L 
is the Gaussian integer lattice expanded by a factor of uj'. In this case we will 
see in section (4.1.2) that as z travels along the straight line from uj' to uj' + 2uj, 
the point (x, y) — (p(z), p'{z)) moves around the real points of the elliptic curve 

= 4a;'^ — g2X between —K'/ (4\/3) and 0, and as z travels along the straight 
line from to 2w, the point {x,y) = (p(z),p'(z)) travels through all the real 
points of this elliptic curve which are to the right of (^^lO). 

We also note that in the case A 7^ 0, the points of order 3 correspond to 
the inflection points of the elliptic curve. We found that at these points and 
for a positive Cosmological Constant the deceleration parameter tends to — 1, 
leading to the asymptotic value for the Hubble's Constant given by Ac^ = 3i?^ . 
Geometrically, points of order three are the points where the tangent line to the 
cubic has a triple order contact. 



■^A point P is said to have order m if mP = P + P ■ ■ ■ + P = O but m' P ^ O for all 

m summands 

integers 1< m' < m, and O is the identity element of the elliptic curve group. If such an m 
exists, then P has finite order ; otherwise it has infinite order. In order to find points of order 
dividing m, we look for points in the FPP such that mP £ L. (Of course, we send the points 
in L, which are the poles of p, to O). 
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5 Analysis and General Results 



5.1 Introduction 

Before presenting the results for the scale factor R{t) in the general case, in 
which both the discriminant A and the cubic invariants g2 and gs are non-zero, 
we shall present three interesting cases in which one of the above mentioned 
parameters vanishes. Each of these cases has a particular physical interpretation 
and special mathematical properties which, with the exception of A = 0, have 
not been discussed in the literature before. The three cases considered are: 



The presentation and ordering of the cosmological results has been dictated 
by the mathematical properties of the elliptic curve: the mathematics deter- 
mines the nature of the allowed physics. 

5.2 Analysis 

5.2.1 The Case of Vanishing Discriminant, A ~ 

There are three limiting cases in which the discriminant A of the cubic equation 
vanishes. This occurs when at least two of the roots coincide. As one can 
see from the form of the cubic invariants, g2, gs Eg. (^) , the discriminant A 
vanishes when either (or both) of A, M is zero. 

In the third case, A is identically zero for the choice of parameters given by: 



In all three cases the solutions for the cubic equation are no longer described by 
elliptic functions and they lose their modular properties (the WeierstraB func- 
tions degenerate). The resulting solutions are expressed in terms of elementary 
functions and represent the standard textbook solutions of Cosmology p7[ |. 

5.2.2 Class of Solutions with a non-zero Discriminant and One of 

the two Cubic Invariants zero. 

This section examines the two special cases where one of the cubic invariants 
32,53 is zero for a non-zero discriminant. It will be shown that the resulting 



A = 0, 
A > 0, 
A < 0, 



92 ?^ 0, 53 ^ 0, 
92 7^ 0, 53 = 0, 
92 =0, 53 7^ 0. 



(54) 
(55) 
(56) 




Critical ' 



(57) 
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solutions are associated with special elliptic curves and have important mathe- 
matical and physical properties. 

The first case is when (72 = , i.e. K' — (case of a flat universe), 173 = 
— j^A'AP 7^ . Here the discriminant A is negative and the resulting two types 
of solutions are characterised by the sign of A, which in mathematical terms 
correspond to the fixed points in the upper half complex plane (r = ±i + ^z). 

The second case is when 53 — 0, i.e. A' — j^-^ and g2 7^ 0. Here A = 
> and the resulting solutions are characterised by the sign of the curvature 
constant K\ The solutions correspond to the other fixed point in the upper half 
complex plane (r = i) . These special points for r in the upper half complex plane 
have significant symmetry properties which will be discussed in later sections of 
this paper. 



Euclidean Universe (A < 0, g2 — 0, 33 ^ 0) 

For K = 0, A 0, M ^ 0, g2 vanishes while gs and the discriminant A are 
both non-zero. This is the case of the Euclidean Universe where the solutions 
are described by elliptic functions. In this case: 

92^0, 53 = ^A'm2, a = -27532<0, (58) 

the roots are given by: 



ei = — v/"36A'M2, 

62 = ^^-36A'Af2, 

63 = ^^-36A'M2, (59) 
and the cubic equation becomes: 



p'^(z)-4p^(z)-<73. (60) 

For this case the discriminant is always negative and that there are two types 
of solutions depending upon the sign of the Cosmological Constant. Interesting 
physics and mathematics arise when equation ( |60| ) is solved at the zeros ^ of the 
Weierstrafi function, p(uo) = 0, it reduces to: 

*A beautiful formula for determining the zeros of the Weierstrafi function p has been 
obtained by Eichler and Zagier and is given by: p(t, z) = 0'^z = Xt + ii± 

U + i^^^ + 144.«v/6 Tit - r)^dt 
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<P ("o) = -53- 



(61) 



For A' > the defining equation becomes: 
p'("o) = iyfal 

leading to the expression for Hubble's Constant given by: 




This is the expression for the Hubble parameter in de-Sitter phase (inflationary 
scenario). At the point of the complex plane that corresponds to zero of the 
Jacobi modular form p, a point of order 3, the flat universe scenario |^ predicts 
the ultimate de-Sitter phase, given by: 



Ac2 = 3F^, (64) 

where Hao is the limiting value of the Hubble's Constant as > oo, and R~ > 
oo. Equation (|6^ ) connects the Cosmological Constant A, with the asymptotic 
expansion rate, itself a dynamical parameter. 

As emphasised earlier, the limiting value (|6j) is achieved at the 3-division 
points of the elliptic curve that correspond to zeros of the p function, while 
the deceleration parameter q given by Eq.([50|) in this limit tends to —1. The 
deceleration parameter is shown in fig.(|^). There q takes values of 0.5 when the 
matter density is maximum, zero at the inflection points and —1 at the points 
where the Cosmological Constant dominates the dynamics of the universe. 

For A' < the defining equation becomes: 



p'M = V93, (65) 
. /A'Af2 M [A^ 

pM - ^V^^TVT' ^^^^ 



^In fact, as we shall see shortly, the asymptotic value Hcx> of the Hubble parameter for 
A > is reached also in models with non-zero K. 
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leading to a Bubble's Constant given by: 



H = ^^\l^- (67) 

In the Euclidean Universe scenario the lattice L is highly symmetrical in 
that it admits the property of complex multiplication (by t — e^"/^) and the 
Weierstraf? function p obeys the functional relationship: 



p(ru) =T-2p(u). (68) 

The above relationship has been verified, numerically,. The associated ellip- 
tic curve has absolute modular invariant j = oQ. 

The scale factor evolution and other cosmological parameters in the special 
case of Euclidean Universe will be discussed in more detail below. 

An interesting scenario, which is also favoured by recent observational data 
(Type la Supernovae, Cosmic Microwave background radiation experiments) is 
the case with A > 0. In the second graph of fig.(^ we have plotted the scale fac- 
tor versus time for the choice of cosmological parameters: A = 10~^^cm~^, M = 
8.5 X 10^'^cm^s-^,K = 0. 



Here we have integrated (^ij) along the diagonal m(a; -I- a;'),TO G R of the 
FPP in the complex plane. Since A < there is one real (negative) root 62 and 
two complex conjugate roots. At the 2-division point (uj -\- lj'), p(w -I- lo') = 62- 
Fig.(||) shows the behaviour of the WeierstraB p function along the diagonal of 
the FPP. As we can see the p{u) function starts from 00, passes through a zero, 
and reaches the negative real root 62 as u varies from 0, through and onto 
w -I- uj' . When u — > vq, R{t) — > 00, t — > 00, Hubble's Constant goes to the 
limiting value given by, — > this is the inflationary scenario. 



Finally we note that when the Cosmological Constant is negative we obtain 
an oscillatory solution see fig.(0). The scale factor has a maximum value at the 
2-division point lo + uj' m the FPP. At this point the p function is equal to the 
real root , i.e. p{u) + uj') = 62, and the first derivative p' vanishes leading to a 
zero value for the Hubble Constant as expected. 

^The absolute modular invariant is defined as; j = 1728 X % = 1728 X , ■> . We 

discuss in more detail its properties in section 7. 
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Non-Euclidean Universe (A > 0,172 7^ 0,173 = 0) The second special case 
is given by the foUowing conditions: 



K' 

ff2 = ^, 53-0, A-.gf. (69) 

In this case A > and there are three real roots, one of which is zero and the 
other two are are equal in magnitude and opposite in sign. The elliptic curve is 
given by the expression: 



y2 = _ ^^^^ (70) 
or in the Weierstrafi parametrisation: 

p\z)^^p\z)-g^p{z). (71) 

The vanishing of 33 leads to the expression for the cosmological parameters, 
given by: 



""-kip- h""" 

The sign of the Cosmological Constant depends on the sign of the curvature 
coefficient. Also, the elliptic curve has the property of complex multiplication 
(CM), in that the square lattice L admits complex multiplication by i. The 
Weierstrafi function pix) satisfies the functional relationship: 



piix) ^ i ^p(x) = -p{x). (73) 

The absolute modular invariant function for the elliptic curve (|7l| ) has the 
value j = 1728, and since the elliptic curve has the property of CM the absolute 
invariant is an algebraic integer. 

There is a very interesting physical feature: at the point vq that inverts the 
equation p{vo) = — the derivative of the Weierstrafi function satisfies the 
elliptic equation, namely: 



K'Y K'^ f K'\ 8K 



or 



^>°)-^^(-T2 -12 =12^' (7^) 



= (75) 
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at this point the Bubble's Constant becomes: 

n.-{l),'M.-^.-c^^. (76) 

Therefore we arrive at the de-Sitter functional relationship for the Hubble's 
constant for A > 0, and again the limiting value, -ffoo, for the Hubble's Constant 
is recovered Eq. (^4|) . 

In fig.(||) we have plotted the scale factor versus time for the choice of pa- 
rameters: K = 1, M = 3.5 X 10^^cm^s~^, A = jj]^j2 - Here we have integrated 
the Friedmann equation in the region, (lu + lu'^ lu' + 2 lu) of the FPP in the com- 
plex plane, where, p{u) lies between 63 — and 62 — 0. In this "bouncing" 
universe scenario, the scale factor starts at a minimum value, i?min(i), and then 
expands exponentially forever. Interestingly, such models do not suffer from a 
t — singularity. 

The deceleration parameter q, versus time, is also plotted in fig (^. In this 
case the deceleration parameter is always negative see second graph. 



Fig.(|[) (first graph) presents results obtained by integrating eq. (p|) , for the 
same values of the physical parameters A, M, K, in the region of the complex 
plane from to 2uj. In this case a periodic solution is obtained, and in this 
region of the FPP, the Weierstrafi function is positive, p(u) > ei — For 
A < an oscillatory solution is obtained (second plot) of Fig.(^). 

In both cases A = ^Acrit, where Acrit is given by equation ([s^). 

5.3 General Results (A 7^ 0, ^2 7^ 0, 0) 
5.3.1 Approach 

The results will be presented in three different classes determined by the sign 
and magnitude of the Cosmological Constant. For each of the classes arbitrary, 
but realistic, sets of cosmological parameters have been chosen (A, M, K). Here 
the pseudomass and Cosmological Constant parameters always have a non-zero 
value such that g2,93 and A are all non-zero. The types of solutions can be 
classified as: 

Periodic or Oscillating Universe, 

Asymptotic Inflationary Universe, 
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Periodic and Bouncing Universe (for < A < Acrit and K = 1) 

Solving the Friedmann equation exactly leads to a set of interdependent 
cosmological physical quantities: 

Scale Factor 
Age of the Universe 
Hubble's Constant 
Deceleration Parameter. 

For each class of solution we have picked (fixed) values of the cosmological 
parameters that give "good" results for the cosmological values. Qualitative 
similar results to our own for the scale factor have been found in Ref. ||2^ . 

In the following section we present our results for the scale factor, Hubble's 
Constant and deceleration parameter. 

5.3.2 Periodic or Oscillating Universe 

For a negative Cosmological Constant the solutions are always periodic. We 
have chosen a representative set of cosmological parameters given by: 

K^0,±1, Af = 2.0 10^^ cTO^s"^ A = ^10"^^cTO-^ (77) 

The resulting solutions are shown in fig(0). We have repeated the calcula- 
tions for a range of values of the cosmological constant (— 10~^^cm~^ < A < 
10~^^cm~^) and a range of values for the parameter M. The values of the 
parameters in (^^ and in following sections were chosen in order that the re- 
sulting cosmological solutions were as realistic as possible (i.e. not too small 
scale factor, or age of the Universe). 

The solutions show a general trend of increasing scale factor and age of the 
universe as K goes through 1, and -1. As expected, shown in fig.(0), the 
deceleration parameter is always positive. For this choice of the cosmological 
parameters, a "good" set of results, for K = 0, is given by: 

Scale Factor 2,901 Mpc 

Age of Universe 10.01 Billion Years 

Hubbies Constant 49.4 Kms^^Mpc~^ 

Deceleration Parameter 2.26 

5.3.3 Asymptotic Inflationary Universe 

This scenario used the following set of Cosmological Parameters: 

K ^0, ±1, M = 8.5 X 10*^ - 4 X 10''^ cm^s"^ A = lO^^^cm-^. (73) 
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For a positive Cosmological Constant the solutions are asymptotic inflation- 
ary for iiT = 0, —1. For A > Kcrit and if = +1, the scale factor time evolution 
is similar to the cases that correspond to K = 0,-1, where eventually the 
asympotic de-Sitter Universe is reached. 

Interestingly for the case, K = +1 and < A < Kcrit, a periodic or a 
bouncing solution is obtained. This class of solutions includes the special case 
ff3 = 0,(72 7^ discussed earlier in the paper, where A = (see also next 

section). The resulting solutions are shown in fig.(||). Again, the solutions show 
a general trend. Initially matter dominates the evolution of the scale factor, 
it then goes through the point of inflection and finally tends asymptotically to 
infinity. The Hubble's Constant in this last phase is determined entirely by 
the Cosmological Constant Eq.(|^ whose asymptotic value is H"^ — i^^Y'^ ~ 

(^A^^y , for the value of A chosen in Eq.Q. 

Similarly, we observe that the deceleration parameter starts with a positive 
value {q — 0.5), goes through zero at the inflection points and tends asymptot- 
ically to minus one. For this choice of the cosmological parameters, a "good" 
set of results, for if = 0, is given by: 

Scale Factor 3300 Mpc 

Age of Universe 13.6 Billion Years 

Hubbies Constant 66 Km s~^Mpc~^ 

Deceleration Parameter -0.477 



The above cosmological values are in agreement with recently published 
experimental results: scale factor deceleration parameter q ^ and Hubble's 
Constant [|@. 

5.3.4 Periodic (special case) and Bouncing Universe 

This section examines the class of models that lead to the periodic and bouncing 
Universe. These solutions are produced when, K — 1 and < A < Acrit- For 
this case we have chosen a parameter set given by. 



is: = +l, Af = 9 X lO^^cm^s-^^ A = 10"^'^ cm"^ (79) 

The resulting solutions arc shown in the first two graphs of fig.(^. For 
comparison the cases with K = —1,0 and same value for the A, M parameters 
are shown in the last two graphs of fig.(^. 

The bouncing solution is of particular interest in that there is no i? = 
singularity, and the periodic solution, for this special case, is similar to solutions 
with negative Cosmological Constant and positive curvature coefficient. 
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6 Discussion 



A new exact solution of the field equations of General Relativity, in closed 
form in terms of the Weierstraf5 modular form for the Szekeres-Szafron family 
of inhomogeneous space time line element with a non-zero Cosmological Con- 
stant has been derived. Our approach was to integrate the non-linear PDEs 
by bringing them to the form of the defining equation of the elliptic curve by 
appropriate substitutions. By careful choice of parameters and constants the 
homogeneous FLRW limit of the space time line element was recovered. The 
same techniques used for the integration of the field equations resulting from 
the Szekeres-Szafron spacetime were applied for the integration of the Fried- 
mann field equation. Within the homogeneous cosmology paradigm a useful, 
predictive theory of cosmology has been developed. Here knowledge of the cos- 
mological parameters K, A and M leads to the derivation of the geometry of 
the Universe where actual values for the scale factor, Hubble's Constant, decel- 
eration parameter are obtained. 

For the cosmological parameters. A' = 0, A = +lQ^^^cm^^ , M = 8.5 x 
10'*^CTO^s~^, a set of values for the cosmological quantities are shown below: 



Scale factor R{t) 


3300 Mpc 


Age of the Universe 


13.6 Billion years 


Deceleration parameter 


-0.477 


Hubble's Constant 


eeKms-^Mpc^i, 



which are in reasonable agreement with the published experimental values 

However a set of cosmological parameters with a negative Cosmological Con- 
stant also predict a "reasonable" set of values, the essential difference between 
the two is the sign and magnitude of the deceleration parameter (g = —0.48 for 
A > 0, 5 = 2.2 for A < 0). It will be interesting to see how the value for the 
deceleration parameter changes over time and consequently the geometry of the 
Universe, as more experimental data become available. 

The body of the paper shows how the use of WeierstraB functions, and the 
associated theory of elliptic curves and modular forms, provides the mathe- 
matical framework for solving several non-linear equations arising from General 
Relativity. Further we suggest (see appendix D) that many other non-linear 
partial differential equations of mathematical physics could be solved from just 
such an approach. There we give examples from soliton physics. Thus the 
methodology and the underlying philosophy developed in this paper is quite 
generic and suggests a unified approach for investigating and solving various 
physical problems of mathematical physics. 

Weierstrafi functions are important in that they provide the link between 
pure mathematics analysis techniques and the physical problems of theoretical 
cosmology. In particular, for a positive Cosmological Constant A and A' = 
(Euclidean Universe), the discriminant A is negative resulting in one negative 
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real root. In this case the integration path includes a zero point, vq, of the 
Weierstrafi function {p{vo) = 0) which invariably leads to an asymptotically de 
Sitter phase in which the Cosmological Constant dominates the dynamics and 
the fate of the Universe, i.e. ^ = H^.The zeros of the p function in this case 
are points of order 3 (i.e. they are points of the FPP that satisfy 3P = C) 
and belong to the torsion subgroup of the elliptic curve group. The eternally 
accelerating universe (which will exhibit an event horizon) is associated with the 
zeros of a modular form! In the general case of positive Cosmological Constant, 
and for K = —1 and K — +1 {A > Kent), we also reach the dynamical de -Sitter 
phase asymptotically at the points vq of the fundamental domain which solve 
the equation p{vq) = —jj- ^^^^ note that in the special case 53 = 0,(72 7^ 
the Universe recollapses even for a positive Cosmological Constant for K = +1. 

For a negative Cosmological Constant and A < , the Weierstrafi function 
is always positive in the physical region of integration in the FPP, and its values 
lie between the positive real root and +00. In this case the integration path in 
the FPP never includes the points that solve the equation p{vo) = —j2 ^'^^ ^ 
periodic solution is always obtained. The maximum scale factor occurs at the 
2— division points, 62 =p{oJ + uj') of the Weierstrafi function. These points are 
points of order 2 which also belong to the torsion subgroup of the elliptic curve 
group. The value of the positive real root determines the maximum value for 
the scale factor. Similar arguments hold for the positive discriminant case. 

Oscillating universes have a certain poetic almost transcedental appeal and 
many ancient cultures fostered the notion of a universe that periodically died, 
only to rise phoenix-like from the ashes to seed a new creation. As described 
by Herodotus Q: "the sacred bird was indeed a great rarity, even in Egypt, 
only coming there according to the accounts of the people of Heliopolis once 
in five hundred years when the old phoenix dies." For a recent account ofthe 
oscillating case see the book by Barrow [ f7|. 

The theory of Taniyama-Shimura |43, |33[ has interesting ramifications, for 
the special cases of the invariants 32 = 0, 53 7^ and 32 7^ 0,53 — 0, which corre- 
spond to the Euclidean and A=^Acrit solutions. These solutions correspond to 
elliptic curves with the property of complex multiplication and are connected by 
the Taniyama-Shimura theory, to modular elliptic curves. The corresponding 
tori admit extra symmetries and their complex structure belongs to a quadratic 
imaginary field (in fact the complex structure is a fixed point of the modular 
group). Consequently, the solutions are intimatetly connected with the hyper- 
bolic upper half-complex plane. For discussions in the mathematical literature 
of the property of complex multiplication we refer the reader to ^ . 

We believe that the machinery developed in this paper based on exact solu- 
tions of General Relativity should be a useful tool for observational cosmology 
and for investigating in a precise manner the physical implications and mathe- 
matical properties of various theoretical models. Finally, the presently accepted 
model of the Universe corresponds to a solution that admitts complex multi- 
plication ((72 = 0, (73 ^ 0) and therefore has maximal symmetry. We believe 
that this is not coincidential but points to a deep connection between General 
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Relativity and topology. This raises the question of Einstein: "Did God had any 
choice in creating the Universe?" . 
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A Ricci Tensor and Scalar for the Szekeres-Szafron 



Metric 

Below arc the equations for the non-zero components of the Ricci tensor, R^i, 
and R scalar. 



^00 = 2^ + ^ + 2 



dp 
'di 



fdaV 



(80) 



-Rii = 



d^p 



a a 



dp 
'di 

da 
dx 



+ 



+e 



2(/3-a) 



5^ 
dz^ 



+ 2 



dp_ 



dp da 
'di'di 

dp da 
dx dx 



d'^P dp_da 
dy"^ dy dy 



dadp 
dz dz 



(81) 



i?22 = -e 



2/3 



df^ 



dt 



dt dt 



d^P dadp 
dx^ dx dx 



d^P d'^a 
dy^ dy^ 



dp da 
dy dy 



da 
dy 



+e 



2(/3-a) 



d^p — 
dz'^ \dz J dz dz 



(82) 



R. 



'33 



o2a 



d'^a / da 
W \~dt 



^dadp 
'dt'dt 



dx'^ 



da. 
dx 



+e 



+2 



2(c«-/3) 



d^a 
dy"^ 




d^^fdp 
dz"^ \ dz 



dadp 
dz dz 



(83) 
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B Integration of the Non-Linear PDE 

The elliptic integral for time can be defined as: 



t = J[-K{z) + 2M{z)^-^ + ^^^]-^/^d^, 



$-1 



-K{z)^ + 2M{z) + ^^^ 



-1/2 



let us first consider the integral, 



-K{z)^ + 2M{z) + ^^^ 



-1/2 



(85) 



-K{z)^ + 2M{z) + ^^^ 



-1-1/2 



(86) 



Defining a new variable X — enables the integral to be reduced in cubic 
form, given by: 



-2M{z)X^ - K{z)X^ + 



A' 



-1/2 



dX 



(87) 



Now by defining a second variable X = — ^~^m{z)/2'^ ' above integral can 
be written in a Weierstral3 normal form: 



/( 



^ 12 ^ 216 ^ 12 



-1/2 



(88) 



Thus ^ = p(u + e), and e is a constant of integration. The cubic invariants 



are: 
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Finally, substituting back into our previous equations, gives: 
and 



t = I Mu = 



p{u + e) + Kiz)/12 
M{z)/2 



du. 



(89) 



(90) 



(91) 



p{u + e) + ^ 

C Derivation for the Energy Density expression 

In this appendix we derive the expression for the energy density using the first 
field equation given by, 



-2 (a + 2/3 + + 2/3 



~«(P + 3P), 



(92) 



Ltq — ZLt^ — Lrg 



$($' + 

2 



($' + $1/') 



($' + $i^') 



Expression 



= 2e'^[eP^]- {^\,, + v'<^,t + 2^{'^' + <^v')} (93) 
can be written as follows: 



$;,, + i.'$,,t + 2^($' + a>z.') 



2M{z] 

$3 
1 K 



kP 
— 



2 



-2M{z) 1 K 



P 



rfi - kP 



1 KV 



2 $2 

3kP 



(94) 
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Eq.(p^) can be written using (p^), 



-k{p + 3P) 



Finally we get, 



(96) 

D Differential equation of the p function, and 
non-linear equations of physics 

The Weierstrafi p{u) = p(u, r) = p(u|w,a;') is an even meromorphic elliptic 
function of periods 2a;, 2uj' (i.e. p(u + 2uj) = p{u) = p{u + 2uj') for all complex 
numbers u), which is of order 2, has a double pole at w = and the principal 
part of the function at this pole being u"^. Thus the Weierstrafi p function, for 
T G 3 (upper half-complex plane) and u G C is defined as follows: 

The Weierstrafi function p (as it can be proven from the definition of the func- 
tion) has non trivial modular transformation properties. It is a meromorphic 
Jacobi modular form of weight 2 and index and signature 1 1|2^ , namely: 



+ = {cT + drp{u,T), ( I 5 ) e5i(2,Z),(98) 



CT + d CT + d) J i \ ^ y c d 

p{u + \T + ^Ji,T) = p(u,t) ((AM)e^')- (99) 

Next we note that the Weierstrafi elliptic function p obeys the following 
important differential equation: 

(p')2=4p3-52(T)p-53(T), (100) 



By differentiating (100) twice we obtain 

p"'(z) = 12p(z)p'(z). (101) 
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The first equation (IOC), is the equation of an clhptic curve, a Riemann 
surface of genus 1. The second differential equation that the p function satisfies 
(101) is the famous Kordeweg-de Vries (KdV) non- linear partial differential 
equation of soliton physics (in a time independent form): 



Ut = Uxxx — ^'2uUx, u = u{x, t) (102) 

In fact the techniques developed in this paper for solving the gravitational 
field equations may have applications to other non-linear differential equations 
of mathematical physics. In particular non- linear PDEs of solid state physics be- 
sides (KdV), like its two dimensional generalization, the Kadomtsev-Petviashvili 
(KP) iQ, the Sine-Gordon as well as the non-linear Schrodinger equation of 
Quantum mechanics can readily be solved by elliptic theta functions associated 
with Riemann surfaces |Q . The key property is that the theta functions asso- 
ciated with a period r (or the matrix for a higher genus Riemann surface) 
that comes from a Riemann surface posses special function theoretic properties. 
One of the most striking properties is that these special 'i?'s satisfy simple non- 
linear partial differential equations of fairly low degree, similar to the differential 
equation that the Weierstrass, p, function obeys. All the above interesting con- 
nections have led us to propose a conjecture, namely: That all non-linear exact 
solutions of the resulting non-linear partial differential equations (PDEs) of Gen- 
eral Relativity with a non-zero A are described by theta functions '0{z\T),i!){z\n) 
associated with Riemann surfaces. Our exact solution given in terms of the 
Weierstral5 p function, of the fairly general class of inhomogeneous Szekeres- 
Szafron metrics already provides strong evidence. Other non-linear equations 
of General Relativity whose solutions described by WeierstraB functions have 
been obtained by (Lemaitre Omer |Q, Zecca [^). 
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Figure 3: Deceleration parameter and scale factor in the Euclidean Universe, 
with parameters A = IQ-^'^cm-^, M ^^.5 x IQ'" cm^ s''^ , K = 0. 



10 20 " — 3Vr 40 50 60 



iii(a)+(y'),iiieQ 



Figure 4: Values of the Weierstra fi pfunction along the diagonal {lv + lo') of the 



Bouncing Universe, K= 1 ,A=(1 /1 8)K'M"',M = 1 0'^cm's"' 
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Figure 5: Scale factor and deceleration parameter versus time in the bouncing 
universe for K = +1,M = 10^^cm^gg^,A = j^K'^M^, in the special case 
53 = 0. 
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Figure 6: Periodic solution for A > (first graph) and A < (second graph) in 
the special case = 0. 
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Figure 7: Periodic solutions for A < 0. 
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Figure 8: Asymptotic inflationary cosniological scenarios for positive Cosmofog- 
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Figure 9: Scale factor vs time for < A < Acru- 
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